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A first principle explanation of the origin of intermittency and nonlinear struc- 
ture formation in the Lagrangian velocity increments of a turbulent flow is presented 
in the context of a scale invariant analytical formalism that is being developed re- 
' cently. The copious generation of power laws and nonlinear exponents in the struc- 

ture functions are shown to follow quite naturally in the present formalism. 
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^ '• 1 Introduction 

o ■ 

-— i \ Understanding the emergence and subsequent control of turbulent flows in fluid motion 
described by the Navier-Stokes equation is of considerable interest recently. Starting from 
the seminal contributions of Richardson and Kolmogorov, there have been great advances 

^ ■ in various directions [I], for instance, the statistical explanations of energy casceding 
c5 ' and dissipation, emergence of intermittency pioneered by Kolmogorov [3] and Obkubov 
citeobh vis a vis the analysis of dynamic equations of hydrodynamical forces given by 
the Navier Stokes equation based on Shell models [2], Stochastic Langevin equations [5 J 
and so on. An interesting area of investigation is the exploration of the relationship of 
turbulence and chaos in deterministic systems. Although the question of existence of 
chaos in Navier Stokes equation is still open, there have been some recent advances in the 
existence proofs of simpler equations such the Euler equation [6]. 

Origin of intermittency in the dynamics of turbulence is a very active field of current 
research [5j [7J [8j [9]. Such studies are recently being addressed mainly mainly in the 
framework of the Lagrangian approach: experimentally, one follows the motion of a single 
tracer particle and the relevant quantity of study is the increments in time of its velocity 
fluctuations A T v(t) = v(t + r) — v(t), where v(t) = v(x ,t) parametrises the velocity of 
the fluid particle initially at Xq along a Lagrangian trajectory j5]. It is suggested that 
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the anomalous scaling in the Lagrangian velocity increments could be traced back to 
the existence of long time correlations in the particle accelearations that relate to the 
hydrodynamic forces driving the particle motion in the turbulent flow. These long time 
correlations are expected to be of fundamental importance in understanding both the 
dynamical origin of turbulence and the associated statistical modelling. 

A key issue in the intermittency research, in the context of a turbulent flow, is to look 
for a dynamical explanation of the anomalous scaling exponents in the structure functions 
of the velocity increments < |A T v(t)| <? >~ because of the extreme non-gaussianity 
of the corresponding probability distributions. A linear exponent of the form £(g) = q/2 
corresponds to the Browninan motion and nonintermittent Kolmogorov flow. A nonlin- 
ear behaviour of £(g) reveals the presence of multifractality leading to intermittency [7]. 
Thus intermittency is linked to the nature and distribution of small scale structure forma- 
tions on the fluid particle trajectory leading to intense bursty activity over a more quiet 
background. 

In this paper, we present a novel scale invariant dynamical (more generally, analyti- 
cal) framework which can support such nonlinear structure formation leading to localised 
bursty activity in a very natural way pm EH EEl [131 03] ■ We also identify the scaling 
exponents with the Hausdorff dimensions of a class of Cantor like fractal sets those are 
shown to arise spontaneously on the trajectory of a Lagrangian particle along the turbu- 
lent flow. This approach is likely to offer a general framework for addressing key issues 
related to chaotic and turbulent motions in geometrically and intuitively much simpler 
manner compared to those available in the literature. It should also offer new insights into 
the emergence and evolution of multifractal structures in various fields such as turbulence, 
finance, web traffic and so on. 

To highlight the actual status of the present work against the current activities [5j [7J 
[8], [9] in the areas of intermittency and multifractality in fully developed turbulence, let 
us mention that understanding the emergence and dynamics of intermittent multifrac- 
tal structures in turbulent flows are considered to be an open issue. Current researches 
mainly aim at developing analytic concepts under various special assumptions [7J [8] and 
computational techniques [5] in analysing relevant questions and interpreting experimen- 
tal facts. This work now presents a novel analytic mechanism for generating copiously 
intermittent and multifractal structures by exploiting a new class of nonsmooth solutions 
to linear differential equations (Sec. 2) that becomes available when the possibility of a 
new level of inversion induced nonlinear increments for a dynamical variable is formulated 
rigorously (Sec. 3). The relevant concepts and results are formulated and developed in the 
context of a Langevin stochastic differential equation when the noise term is assumed to 
have a special scale invariant form (Sec. 2). The random fluctuations in the stochastic 
noise term is then related to the inversion mediated jump distributions on a Cantor like 
fractal set leading to nontrivial velocity increments. A one parameter class of such Cantor 
sets is shown to arise naturally in the present scale invariant formalism applied to the tur- 
bulent flow. The bursty intermittent behaviour is then interpreted as due to occurrence of 
Devil's staircase functions as ultrametric weights (valuations) assigned naturally to such 
multifractal Cantor sets. As a consequence, this extended scale invariant analysis and 
the associated dynamical framework turns out to offer a natural framework to study dy- 



2 



namics of trubulent flows (Sec. 4) and also of other complex systems. In short, this paper 
offers a new approach towards ushering in a paradigm shift in addressing key questions 
related to turbulent flows. The standard classical analytic framework might very well be 
extended to a scale invariant formalism that has been shown here to offer more natural 
interpretations to observed facts such as intermittent structure formations in turbulence. 
The present work is, however, mainly of theoretical interest and does not yet make any 
contact to experimental data. 



2 Motivation 

The formulation of the scale invariant analysis is motivated by the possible existence of 
a class of nonsmooth solutions of the simplest ode of the form 



dv 1 

Tt = ~¥ V (1) 

The formalism of the scale invariant analysis is developed originally in the context of 
an extended real analysis incorporating dynamically active infinitesimals living in a non- 
archimedean space jT2]. However, we outline here essential features of the said analysis 
using the language of turbulence. As a consequence, v(t) = v(r ,t) here is assumed to 
represent, for definiteness, the Lagrangian velocity field of a tracer fluid particle in a 
turbulent flow. In ordinary sense, the above equation of course represents a laminar 
phase that relaxes in time T . 

The origin of nonsmooth solutions can now be explained as follows. Let us begin 
by noting that a solution is nonsmooth in the sense that it could at most be realised as 
an approximate solution in the framework of ordinary calculus and also in the classical 
mechanics. However, in an extended non-archimedean, ultrametric space, such a solution 
can indeed be accommodated as a smooth solution. Our aim here is to highlight how such 
nonsmooth solutions could play a significant role in generating intermittent behaviour in 
a turbulent flow. 

Notice that in ordinary calculus, an independent real variable t, considered here as 

time, is supposed to undergo changes (increments) only by linear shifts of the form t — > 

t_ 

t + At, and eqfCQ) is known to have the unique solution v (t) = v(0)e T o . Suppose that a 
localised disturbance is created (by stirring, say) on the initial laminar flow of fluid, say, 
at time T . The subsequent turbulent flow is modelled phenomenologically as a stochastic 
Langevin equation of the form 

— --— v + — (2) 
dt Tq dt 

where forcing noise term is represented as dF = v dl °^ T dt. Such a stochastic modelling is 
expected to describe the fluctuating motion of the tracer particle in a turbulent flow [5] 
with a very low viscosity and hence at a very high Reynold number Re = IqVq/v, where 
Iq is the integral scale, vq stands for the rms (root mean square) velocity fluctuation at 
the integral scale and v is the kinematic viscosity. For a white noise, the above equation 
leads to the Brownian motion. However, in the present approach, the basic Gaussian 
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random variate is assumed instead to be the scale invariant (stochastic) time r, and so 

the noise is log-normal and hence is coloured. The scale invariant time r has a value of 

order 0(1) when the ordinary Newtonian time t is of 0(1). However, r must approach 

to infinity as t — > oo. The above ansatz for the noise term now could be interpreted as 

revealing the existence of long time correlations in the time variable t itself due to the 

induced fluctuating variable r, thus having the novel representation t ~ T (l + 77 + log r), 

rj being a growing (0(1)), ordinary (real) variable, and T — >■ 00. Notice that t reduces to 

the ordinary time variable viz, t ~ Tq(1 + rj) when To ~ 0(1). This long time correlated 

fluctuation then, in turn, will be inherited by the underlying velocity fluctuations as well. 

As a consequence, noise term is incorporated into the very definition (ie, the intrinsic 

structure) of time itself (cf. next section for more precise justifications). This is one of the 

t_ 

salient features of the formalism. Removing the laminar part v = v(0)e T o , t ~ T (l+r)), 
the renormalised velocity fluctuation v = vv^ 1 is now given formally by the scale invariant 
equation 



Eq©, although should be interpreted as a (coloured) stochastic equation, however, in 
the following is treated as an ODE (deterministic system) written instead over a positive 
measure Cantor set : a fattened real number system 1Z, each point of which is replaced 
by a zero measure infinitesimal Cantor set C p belonging to a one parameter p family of 
such sets. The elements of the Cantor set C p o, that replaces the ordinary of R, are 
new infinitesimally small elements in 1Z, known as soft zeros having dynamic properties. 
The scale invariant ODE ([3]) is well defined over the gaps, of all possible sizes, of such 
a Cantor set [12]. The randomness in the fluctuations is therefore transferred to the 
dynamical formation of Cantor like sets in the back ground of a turbulent flow. 

Now to formulate motion over a Cantor set rigorously, we make an assumption that 
the scale invariant time r undergos increments by random inversions of the form r — > r _1 , 
in the neighbourhood of r = 1 (say) . Let r_ = 1 — rj(t) and r + = 1 + fj(t) such that the 
inversion rule t_t + = 1, or more generally, r + ri +a = 1, a being a real parameter indicating 
the size of the increment (jump), is satisfied. The inversions are realised randomly in the 
sense that a fluid particle at a site r_ jumps over instantneously to the site r + when r_ 
and T + are separated by a gap of size a. Because of small scale (Cantor set like) structures, 
the exact specification of r + in the fattened real number set TZ is also not precisely known. 
The inversion rule then defines a small scale multiplicative structure close to r = 1, that 
arises dynamically as the ordinary time t approaches asymptotically the limiting value 
00 (cf, Sec 3 for details). The multiplicative scaling exponent a in r + = rZ x ~ a is an 
ultrametric valuation defined by 



line segment neighbourhood of T is now extended to a positive measure Cantor like 
set C p , realised as an ultrametric space with the absolute value a. Transition between 



dv = vdlogr 



(3) 





as t — > 00. As a consequence, and also because of scale invariance, the ordinary connected 
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granular like points of such a Cantor set, realised as end points of the gap, are mediated 
by seemingly stochastic inversion induced jumps (of sizes determined by the associated 
gap sizes) as defined above, so that the corresponding multiplicative jump differential 
measure is defined by Ajt = loga(t), t — > oo. Notice further that the valuation a is a 
locally constant function defined by da{t) = which is a direct consequence of the fact 
that d\og\ogr + = d log log tT 1 , so that a(t) has a globally intermittent behaviour given 
by a Devil's stairecase function (i.e. a Cantor function) indicating significant fluctuations 
in the neighbourhood of a point of the ultrametric Cantor set (ie. the set over which the 
valuation undergoes actual variations). However, a(t), although it is constant over the 
scale of t, ie, ^ = 0, its intrinsic variability is revealed in the smaller log log t scale, ie, 

dr 

^gt- 1 -^— = r. (5) 
alogt 1 

To summarise the main advantage of our analytic approach over others is the dynamic 
generation of infinite numbers of cascaded smaller scales of the form < l n < X 2 j^ for 
every pair of bounded intervals of the form < \\ < A 2 << 1 in a scale invariant manner 
and thereby assuring formation of a host of Cantor like sets, all of which are distributed 
following certain multifractal measures. 



3 Formalism 

3.1 Revisiting Ordinary Calculus 

Let us recall that analysis of motion, from simple and/or laminar motion to more complex 
chaotic and/or turbulent flows, ultimately rests on the basic concepts of differential cal- 
culus. Let us begin first by showing how a slight variation of conventional treatment and 

meaning of the statement lim t = 0, t being a real variable can lead to potentially new 

t— >o 

results which would be significant for justifying the subsequent scale invariant formalism. 
The ordinary meaning of the above limiting statement is, of course, that the real variable 
t approaches linearly and uniformly with constant rate 1 to the limiting value remaining 
always on the connected line segment [0,t]. 

This continuous and connected flow may, however, be broken by introducing a scaling 
of the form f(t) = at, < a < 1 and subsequently iterating infinite number of times. 
The original length (Lebesgue measure) t now reduces to successively by jumps as a n t, 
as n — )> oo. Conventionally, one feels that these are only two possibilities with unique end 
result. 

Surprisingly, however, there exists still another nontrivial possibility. To every t and 
a satisfying above constraints, one is free to choose an arbitrary finite p : 1 < p < a^ 1 
such that the continuous linear measure of the flow t — > is preserved at every level n, 
following a measure conservation principle 

t = {ta nlo ^ 1 }xp nl ^ T (6) 

The scale invariance of the above equality is obvious. The linear flow (in the lhs), as 
it were, is transformed into a scale invariant nonlinear flow in logarithmic scales. As a 
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consequence, the reduction of the linear measure 1, at each scale by the factor 

a n\ og t ig 

corrected multiplicatively by creating a space in a verticle direction containing the fattened 
variable T > 1 in the logarithmic scale and as a consequence, in the limit n — > oo, 
and hence as t — > 0, the linear variable t is transformed into a nonlinear flow pattern 
characterised by the scaling law T _1 = t~ , where s = j^f^r is the fractal (Hausdroff) 
dimension of a Cantor like set that may be assumed to have been formed spontaneously 
in the neighbourhood of 0, because of the limit process. The precise dynamical setting of 
realising a cascade of scales of the form a n and the associated (daughter) multiplicative 
scales p n will be presented in the next section. We note that, because of the inversion 
rule, the original variable t on its approach towards 0, would get reversed (ie, experience 
a bounce) close to 0, so that the corresponding scale invariant infinitely large fattened 
variable T s = T -1 = tT^ 1 will live and vary on a Cantor set of fractal dimension s = 
1Pg iogp " This, in turn, would have a very important consequence: the value is actually 
never reached; the linear flow along the positive t axis, in the interval (0,t), t I, will be 
transferred in a scale invariant manner, to a new branch in the inverted variable T s ~ t s , 
by inversion induced jumps. The subsequent (linear) flow of T s toward will then again 
experience a bounce, and so on leading to a very intricate multifractal tree like structure. 
The multifractality follows from the existence of a one parameter family of multiplicative 
scale factors p(a) = pa depending on a. 

To have a more concrete appreciation of the above analysis, let us suppose that points 
of the line segments [0,1] have the structure of tiny infinitesimally small beads aligned 
horizontally. Under an application of the contraction the length of the original segment is 
reduced to a. How do the infinitesimal beads align on the contracted segment? Since no 
part of the original segment is removed, the beads would climb one above other to create 
a space dynamically; and this process would continue in the limit n — > oo, and ultimately 
one realises a fractured fractal like structure in the neighbourhood of 0, but in the vertical 
direction. 

In the following section, we briefly present the formal structure of the scale invariant 
analysis accommodating the above novel structure in a rigorous manner. 

3.2 Scale Invariant Analysis 

We state the basic definitions introducing scale invariant dynamic infinitesimals, and 
the associated non-archimedean norms and also outline briefly the natural emergence of 
intermittent structures induced by Cantor's stairecase functions (for details, refer to [T2] ) 
and their possible role in the asymptotic behaviours of a dynamic variable related to a 
turbulent flow. 

By dynamic infinitesimals we mean a class of quantities living in a well defined space, 
which are arbitrarily small in comaprison to any finite element of the space, however, might 
have nontrivial 0(1) effect on the long time evolutionary properties of such elements. For 
definiteness, the underlying space is assumed to be the set of real numbers R, and any 
finite real variable is denoted by x G R. In a dynamical problem the variable x is replaced 
by the ordinary time t. 

Definition 1: Let x £ / = [0, 1] C R and x be arbitrarily small, i.e., x ^ 0, but, 
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nevertheless, x — > + . Then there exists 5 > and a set of (positive) relative infinitesimals 
x relative to the scale 5 satisfying < x < 5 < x and the inversion rule x/5 oc 5/x. The 
associated scale invariant infinitesimals are defined by X = lim x/S. 

As an example, for a continuous variable x approaching + and considered itself as 
a scale, a class of the relative infinitesimals are represented as a; oc x l+l (l + o(x)), < 
/ < 1, corresponding to the real variable x l ~ l , so that the corresponding scale invariant 
infinitesimals are defined by the formula X = \x l +o(x m ), m> I. Choosing the parameter 
A from an open set U C (0, 1), we get a class of relative infinitesimals x belonging to an 
open subset of (0,5), all of which are related by the inversion rule to the real number 
x i-i f Qr gj ven x an( j i Moreover, allowing A to vary from a countable number of disjoint 
open intervals U r , r = 1, 2, . . ., one can introduce a countable, disjoint class I r of relative 
infinitesimals, all of which are related by the inversion rules x r /5 = X r 5/x, where A r G U r . 
Notice that \JI r C (0, 5). It follows therefore that the infinitesimals actually reside in the 

r 

gaps of an arbirarily chosen Cantor subset of the open interval (0,5). As a consequence, 
the nonlinearity induced via the inversion rule in the apparently linear limiting motion of 
x naturally leads to a large family of Cantor like fractal sets. 

Notice that in the limit 5 — » 0, the set of relative infinitesimals apparently reduces to 
the singleton {0}. However, the set of scale invariant infinitesimals may nevertheless be 
nontrivial. The set of relative infinitesimals, denoted 0, has the asymptotic representation 
= {0,5X r }, as 5 — > + . For definiteness, the ordinary zero (0) is called the stiff zero, 
when the nontrivial infinitesimals are called soft or dynamic zeros. The ordinary real line 
R is then extended over R = {x : x = x + 0, x G R}, which as a field extension, and 
because of the Frobenius theorem, must be an infinite dimensional nonarchimedean space. 

Definition 2: Non-archimedean Norm The mapping v : — > I + = [0, 1] 

<f>(x) = \\x\ \ := lim log*-i X~\ x = 5X G (7) 

5->0+ 

together with 0(0) = defines a non-archimedean (ultrametric ) norm on the set of 
infinitesimals 0. Infinitesimals weighted with above absolute value are called valued in- 
finitesimals. Clearly, the ordinary Euclidean norm for a dynamic infinitesimal would 
vanish trivially because \x r \ = 5\X r \ = 0, in the limit 5 — > 0. 

Let us recall that an absolute value is non-archimedean if the triangle inequality is 
replaced by the stronger ultrametric inequality v(a+b) < max{t> (a), v(b)}. Ultrametricity 
now tells that the dynamic infinitesimals are indeed awardednontrivial values of the form 
4>(x r ) = a r p~ rs for p = p(5) > 1 and a r > being an ascending sequence continuously on 
the closure of (J I r , so that the valuation (f>(x) has the structure of an intermittent Cantor 

r 

function associated to a Cantor set with fractal dimension s that has formed dynamically 
in the infinitesimal neighbourhood of 0. Since the ultrametric parameter p generates a one 
parameter family of scale factors corresponding to a family of Cantor sets, the nonlinear 
structure that emerges out from the above analysis of course corresponds to multifractality. 

One of the most important features of dynamic infinitesimals (by inversion, infinitely 
large scales) is that these can actually lead to significant influences in the long time asymp- 
totic properties of an evolving dynamical system. As a result, a laminar flow can indeed 
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be transformed into a complex turbulent flow over a sufficiently long time scale and also 
possibly in the presence of a slight disturbance. 

To see the precise manner how the small scale nonlinear structure inhabited by the 
dynamic infinitesimals can induce an 0(1) variation in the long time properties of any 
system, let us first note that the dynamic infinitesimals introduces a sort of a deformation 
into the standard structure of finite real numbers x through the nontrivial valuation in 
the form 



where the valuation 0, realised as a Cantor's stairecase function, satisfies eqflH]). The scale 
invariant factor X± = X±/x corresponding to the fattened variable X, by assumption, 
undergo transition by inversions, so that the scale invariant increments in X may be 
defined as 



which matches with the linear increment 5iX = \X — 1| = 0(x) logx -1 when x — > 0. The 
corresponding linear increment for the fattened variable has the form SiX = \X — x\ = 

0(5;):cloga;- 1 as x — > 0. Notice that these are new level of increments which arise out of 
inversion mediated variations, that in the present formulation are activated in presence 
of dynamic infinitesimals leading to small scale structure formation which in turn would 
create long time correlations and very large intermittent fluctuations [5] . 

4 Applications to Turbulence 

To apply the above analytical results to a turbulent flow, let us notice that the stan- 
dard multifractality assumption suggests formation of cascades of eddies between the 
integral and dissipation scales converging towards a family of fractal sets Ch with fractal 
Hausdorff dimension D(h) for each scaling exponent h from a bounded interval of the 
form (h m in, h max ). The corresponding increments of the Lagrangian velocity fluctuations 
5 T v = v(t + r) — v(t) are assumed to follow the scaling law 5 T v ~ v (r/T Q ) h for each 
scaling exponent h from the bounded interval. Here, vq and r correspond to the integral 
scale characteristic rms velocity fluctuation and the corresponding integral time scale re- 
spectively. Notice that the derivation of such power laws is impossible in the ordinary 
(classical) framework of Navier-Stokes equation. Clearly, the scale invariant analytical 
formalism now offers us with a mechanism to realise this multifractal scaling dynamically 
in a most natural manner. 

Notice that one can identify the scale invariant velocity increment 5 t v/vq in the tur- 
bulent flow with the increment of the corresponding scale invariant fattened velocity S S V 
following the incremental law fl9]) in the limit of sufficiently small (stochastic) time in- 
crements < r « t , so that one obtains 5 T V ~ 0(f) log r -1 . Now, recalling that 
the valuation enjoys the structure of a Cantor Stairecase function, it follows that 
(S T v)/v := 5 T V ~ 0(f) log r^ 1 scales as (S T v)/v ~ (t/to)' 1 ^ f° r an e > (that arises 
form the logarithmic factor), and h being the Hausdorff dimension of the Cantor set Ch 



x — > X± = x ■ x 
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(9) 
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that emerges dynamically in association with the dynamic infinitesimals and the scale 
invariant limiting processes in the present formalism. The above scaling also reveals that 
the velocity undergoes variations only over such family of Cantor sets. Let us emphasize 
that this power law velocity increment follows not from any hypothesis , but presented here 
as first principle derivation from the scale invariant approach. Notice that the set Ch is 
the set that arises in association to the velocity field fluctuations, when Ch is the set that 
is realised in physical space of cascaded eddies in the fluid medium. 

Adopting the well accepted analysis available in literature [U El [7], one then immedi- 
ately concludes that the structure functions of the velocity fluctuations S q (r) = E(\5v\ q ) ~ 
(t/to)^ q \ where the scaling exponent is a nonlinear concave function of q, and can be 
defined following Ref.pQ as a Legendre transform of D(h): £(q) = inf [qh + 1 — D{h)\. 

h 

Here, D(h) > is the fractal dimension of the set Ch on which cascades of eddies accu- 
mulate when the indefinite interations are realised internally as t — > oo. Notice, also that 
1 — D[h) actually stands for the codimension of the limit set Ch in the 1 dimensional time 
axis. The derivation of the scaling exponent £(g) in the context of multifractal random 
walk models [5j [7J will be considered elsewhere. 

To summarise, we have outlined a first principle derivation of nonlinear multifractal 
scaling exponents £(q) of Lagrangian velocity fluctuations. The analysis is based on a 
simple reformulation of the definition of a limit process exposing how small scale structures 
are realised in the neighbourhood of a disturbance. When the laminar flow is stirred 
vigorously, the velocity v(t) = v(ro,t) of a fluid particle at ro at time t would have a 
spectrum of fluctuating values v (t + r) for scale invariant < r << 1 living in a Cantor 
like set Ch with fractal dimension h. Instead of relaxing quickly to the laminar phase, these 
fluctuations will exist as stable structures corresponding to the scaling exponents £(q). 
Conversely, the energy injected at the integral scale by stirring cascades through smaller 
scales following nonlinear convective terms of the Navier Stokes equation. These spatial 
fluctuations then induces Lagrangian velocity fluctuations in the time domain which can 
be interpreted (following the present analysis) as generated from intrinsic, scale invariant 
stochastic variations in the time t itself. As a consequence, multifractal eddy formations 
can be considered as physical menifestations of dynamical multifractals those are shown 
to exist in the infinitesimal time scales. More detailed analysis and interpreting £(q) in 
the light of experimental data will be considered separately. 

Finally, we note that this scale invariant analytical formalism might also be crucial in 
elucidating the nature of the plasma turbulence observed in the space plasma and high 
temperature fusion plasma devices where the turbulence is anomalous. A detailed analysis 
of this aspect will be reported in a separate publication. 
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